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ABSTRACT 



Context. 

Aims. We create a catalogue of simulated fossil groups and study their properties, in particular the merging histories of their first- 
ranked galaxies. We compare the simulated fossil group properties with those of both simulated non-fossil and obserx'ed fossil groups. 
Methods. Using simulations and a mock galaxy catalogue, we searched for massive (> 5 x 10'^ /r'Alo) fossil groups in the 
Millennium Simulation Galaxy Catalogue. In addition, attempted to identify observed fossil groups in the Sloan Digital Sky Survey 
Data Release 6 using identical selection criteria. 

Results. Our predictions on the basis of the simulation data are: (a) fossil groups comprise about 5.5% of the total population of 
groups/clusters with masses larger than 5 x 10'^^ /r' Alo- This fraction is consistent with the fraction of fossil groups identified in the 
SDSS, after all observational biases have been taken into account; (b) about 88% of the dominant central objects in fossil groups are 

elliptical galaxies that have a median R-band absolute magnitude of 23.5 - 5 log h, which is typical of the observed fossil groups 

known in the literature; (c) first-ranked galaxies of systems with A1 > 5 x 10''' h^^Mo, regardless of whether they are either fossil 
or non-fossil, are mainly formed by gas-poor mergers; (d) although fossil groups, in general, assembled most of their virial masses 
at higher redshifts in comparison with non-fossil groups, first-ranked galaxies in fossil groups merged later, i.e. at lower redshifts, 
compared with their non-fossil-group counterparts. 

Conclusions. We therefore expect to observe a number of luminous galaxies in the centres of fossil groups that show signs of a recent 
major merger 

Key words. Methods :N-body simulations- Methods: statistical- Galaxies:clusters:general- Galaxies:evolution 



1. Introduction 



I Jones et all (l2003h identified fossil groups as spatially extended 
X-ray sources with an X-ray luminosity Lx > lO'*" h^^ erg 5"' 
whose optical counterpart was a bound system of galax- 
ies with AM 12 > 2 mag, where AMn was the differ- 
ence in absolute magnitude in R-band between the bright- 
est and the second brightest galaxies in the system within 
half the projected virial radius (r,,,>)- The dynamical masses 
of the systems studied so far are comparable to those of 
rich c l usters (~ 10'^ - 10'^ h'^JVL) dMendes de Oliveira et al] 
20061: ICvpriano et all 120061 : iMendes de Oliveira & et af] 120081: 



Khosroshahi et alj 20061) . Fossil groups may be of considerable 



importance as the place of formation of a significant fraction 
of all giant ellipticals. Beside minor differences in the defini- 
tion of fossil groups, their incidence rate was estimated by ob- 
servational, ana l ytical , num erical, and s e mi-an alytical analyses. 
IViMilinin et"!!] Il999h and iJones et al.l ( l2003b stated that fos- 
sil groups re present (8-20)% of observe d systems in the same 
mass range. Ivan den Bosch et al] (l2007h used the 2dF Galaxy 
Redshift Survey data to measure a fossil frac tion of 6.5% among 
roups with masses (lO'^ - lO'"*) /z 'Alm. [Milosavhevic et al] 
20061) estimated analytically that fossil groups represent 5-40% 
of groups with masses in the range ~ 10'^ - lO''* Ii^^Mq, while 
the percentage decreased to 1-3% for groups of mass larger than 
10'"* /i^'AIq, the latter result having been confirmed using the 



photometric Sloan Digital Sky Su rvey Data Release 2 (S PSS 
DR2). Numerical simulations by iD'Onghia et alj ( l2005l) pre- 
dicted a higher fraction of fossil systems (33%) amongs t group s 
of mass ~ lO'"* /i 'Mq. Ivon Benda-Beckmann et alj feoOTi) . 
also using numerical simulations, found that 24% of groups 
with masses in the range (1 - 5) x 1 0'^ h^^Aig were fossil 
groups. From semi-analytical models, ISales et akl ( |2007|) esti- 
mated that fossil groups represent (8 - 10)% of groups with 
masses (10'^^ - 10'^) h-^Me, while iDariush et al] jmm stated 
that ~ 13% of groups in that mass range were fossil systems and 
predicted that this percentage decreased to 3-4 % for X-ray rich 
systems. 

A natural question is whether the large magnitude differ- 
ence between the first and second ranked galaxies (AMn > 2 
in the R-band), characteristic of these groups, implies that they 
are a distinct class of objects or if they instead represent a tail 
of the cluster (hereafter non-fossil group) distribution. To inves- 
tigate this question, D'Onghia et al. (2005) used high-resolution 
N-body/hydrodynamical simulations to compare four simulated 
fossil and eight non-fossil groups, all of virial masses close to 
1 X 10'"* /i"' Mq. They found that the values of the magnitude 
gap, AMi2, for the 12 systems, were correlated with the halo as- 
sembly time, suc h that fossil groups as sembled earlier than non- 
fossils. Similarly. iDariush et al.l ( |2007|) concluded, by the study 
of fossil groups in the Millennium Simulation, that fossils as- 
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semble a higher fraction of their masses at higher redshifts than 
non-fossil groups. The most accepted scenario for fossil groups 
is then that they are not a distinct class but, are instead, examples 
of groups/clusters that collapsed early. 

Although fossil groups in general do not appear differ from 
galaxy clusters of similar masses, except for their earlier times 
of formation, it was realised, observationally, that the first- 
ranked galaxies in fossil and non-fossil groups differ in some re- 
spects. First, their shapes differ; while fossil, first-ranked, galax- 
ies are often disky, brig htest cluster galaxies are often boxy 
dKhosroshahi et al]|2006l) . Second, their stellar populations are 
dissimilar: fossil group first-ranked galaxies are n ot as old as 
brightest cluster galaxies dde la Rosa & et aL 20081). The s e find- 
ings motivated us to revisit the study of iDariush etaP (l2007h 
of fossil groups in the Millennium Simulation, but now fo- 
cusing on the properties of the first-ranked galaxies (as op- 
posed to the complete system). By searching in the Millennium 
Simulation, we created two samples, one of fossil groups with 
AM12 > 2 (in the R-band) and a second control sample of non- 
fossil groups with systems with AM12 < 0.5. Both the sim- 
ulations themselves and a mock catalogue were used to com- 
plete these searches. In addition, we searched for fossil groups 
in the Sloan Digital Sky Survey D ata Release 6 (SDSS DR6) 
(lAdelman-McCarthv & et"aDl2008h using the same criteria, to 
compare with the mock catalogue as a test of the semi-analytic 
model and also to examine observationally the results of our 
searching algorithm. 

The layout of this paper is as follows. In Sect. |2]we briefly 
describe the galaxies in the Millennium Simulation, and the 
search for simulated fossil and non-fossil groups. Section|3]con- 
tains a comparison between fossil and non-fossil groups, and 
simulated and observed fossils. In particular, we discuss the im- 
plications of these results for the evolution of first-ranked galax- 
ies in fossil groups. In Sect. |4] we describe the construction of 
a mock catalogue and the procedures of group identification. 
We also perform an identification of fossil groups in SDSS DR6 
and compare the results with those obtained from the mock cat- 
alogue. Finally, we summarise the paper in Sect. |5] 



2. Construction of the Simulated Fossil Group 
Sample 

2.1. Dark Matter Particles and Galaxies 

The Millennium Simulation is t he l a rgest c ompleted cosmolog- 
ical Tree- Particle-Mesh (TPM, |x3 (1 19951) ) A?-body simulation 
dSpringel et al. 2005), which evolved 10 bilHon (2160^) dark 
matter particles of mass 8.6 x 10^ h^^Aie within a periodic box 
of 500 /j"' Mpc on a side, using a comoving, softening length of 
5 kpc0 The cosmological parameters of this simulation cor- 
respond to a standard cosmological model (ACDM): Q.,,, - 0.25, 

= 0.75, cTg - 0.9, and li - 0.73. The merging history trees 
were stored for 60 output times separated by time intervals given 
by ln{l + Zn) = n(n + 35)/4200. These are the basic inputs re- 
quired by the semi-a nalytic model. 

We use a run of iDe Lucia & BlaizotI (l2007h 's semi-analytic 
model to extract galaxies with positions, velocities, as well as ab- 
solute magnitudes (in five photometric bands, BVRIK) and stel- 
lar masses, among other quantities. In this model, the branches 
of the halo merger tree are followed forward in time, and sev- 
eral astrophysical processes are included such as gas infall and 



cooling, reionization of the Universe, star formation, black hole 
growth, AGN and supernova feedback, galaxy mergers, and 
spectro-photometric evolution. The final output at z = pro- 
duced ~ 10 X 10^ galaxies with absolute magnitudes Mr - 
5 log It < -17.4 and stellar masses larger than 3 x 10** h^^Mo. 
This version of the semi-analytic model provides an improved 
fit to the bright end of the galaxy luminosity function, compared 
with former model developed bv lCroton et al . (2006). This sam- 
ple is called Millennium Simulation Galaxy Catalogue (MSGC 
hereafter). 

The cosmological parameter set used in the Millennium 
Simulation was that inferred from the first-year WMAP 
(Wilkinson Micro wave Anisotropy Probe Observations, 
ISpergel et al.ll2003l) . Differences in the cosmologi cal parameter 
sets c orresponding to the 1 " and 3'^''-year WMAP dSpergel et al.l 
20_07) (mainly in erg, n) produced a significant delay in 
structure formation in the WMAP3 case, and the number of 
halos drawn for the WMAPl cosmology to be higher. However, 
IWang et akl (|2008|) investigated the implications of this delay 
for the observed properties of galaxies. They compared results 
obtained from a simulation for the cosmological parameters 
of WMAPl plus a semi-analytic model (such as MSGC) 
with those for a simulation for the cosmological parameters 
of WMAP3 plus a semi-analytic model. They found that the 
luminosity functions, correlation functions, and Tully-Fisher 
relations were almost identical in both cases. They stated that 
the galaxy clustering and other observable properties were 
far more sensitive to the galaxy formation physics than to 
the cosmological parameters; the semi-analytic parameters 
then may be able to compensate for the delay in structure 
formation pr oducing galaxy pop ulations at z=0 that agree with 
observations. IWang et aU d2008l) also concluded that substantial 
differences between the models appear at redshifts higher than 
2, which are beyond the scope of this paper 



2.2. FoF galaxy halos 

Since we were interested in studying fossil groups in the MSGC, 
we identified galaxy halos by using a standard method. Groups 
of galaxies in the MSGC were identifi ed by using a Frie nds-of- 
Friends (FoF) algorithm in real space dDavis et alJll985h with a 
linking length of 0.2 of the mean particle density, which corre- 
sponds to an overdensity of 200. We checked that all galaxies 
given in a galaxy halo belonged to the same DM halo. 

Finally, only groups of galaxies with more than 10 members 
were selected We refer to these groups hereafter as FoF galaxy 
groups (or halos). 

For all FoF galaxy groups, we compute the velocity disper- 
sion, virial radii, and virial theorem masses. The value of the 
virial theorem mass was computed to be: 



Mi 



FOF 



(1) 



' The Millennium Simulation, developed by the Virgo Consortium, 
is available at ,http://www.mpa-garching.mpg.de/millennium| 



where Rj, - {^l/Rij^ is the harmonic mean projected sepa- 
ration, given the proiected separations R/j (see eq. [10-23] of 
iBinnev & TremaiMll987h . The virial radius is ri,,> = nR/, and 
the 3-D velocity dispersions were calculated in the MSGC by 
using the peculiar velocities 

2.3. The Sample of Fossil Groups 

To represent fossil groups, we selected FoF galaxy groups with 
masses larger than 5 x 10'"' /i^'AI©, which contained galaxies 
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Table 1. Median Properties of Fossil Groups identified in the 
Millennium Simulation 



Property 


MSGC 


Mock Catalogue 


SDSS DR6 


# 


729 


22 


6 


Mfof [/r'Mo] 


7.5 X 10'^^ 


6.7 X 10" 


8 X 10'^^ 




600 


596 


565 


Tyir [MpC /j"'] 


0.98 


1.0 


1.1 


Mr - 5log{h) 


-23.48 


-23.73 


-22.17* 


M. [/i-'Mo] 


2.8 X 10" 


4.2 X 10" 





Notes: # number of fossil groups, Mfof'- virial mass of the FoF 
galaxy group, 0-30: velocity dispersion of the FoF galaxy group, 
r,.,> ; virial radius of the FoF galaxy group, Mr — 5logh: rest 
frame R-band absolute magnitude of the brightest galaxy within 
0.5r,,„ of the FoF galaxy group, M,: stellar mass of the brightest 
galaxy. 

* In this case it is in the r -band. k+e corrections are computed 
from lBlanton et alj ( 120031) 




10" - 



10" lO's 

Fig. 1. Scatter plot of stellar mass of the brightest galaxy and 
virial mass of the group. Squares correspond to the sample of 
fossil groups, while crosses correspond to non-fossil groups. 



with a magnitude distribution displaying a gap AM12 > 2 (in the 
R-band), when considering objects within a radius of 0.5r,,„ . The 
adopted lower limit to the group mass (5 x 10'^ h^^Mo) ensured 
that the fossil groups chosen were a lso X-Ray fossils, according 
to the work of iDariush et al.1 ( |2007|) . 

In the MSGC, 729 FoF galaxy halos satisfy the fossil criteria, 
which represent 5.5% of the FoF galaxy groups in the studied 
mass range. The median properties of these groups are quoted in 
Table □ 

2.4. The Sample of Non-Fossil Groups 

The principal aim of this pa per is to compare the properties of 
fossil and non-fossil groups. iDariush et al.l ( l2007h analysed the 
differences between fossil and non-fossil systems by concentrat- 
ing on the global properties of the groups. In this work, we intend 
to go further to not only confirm their findings but also extend 
their analyses to the brightest galaxies in these systems. We con- 
structed a control sample of non-fossil groups, using the similar 
criteria to that used for fossil groups in the MSGC but with the 
one difference that the magnitude gap between the brightest and 
second brightest galaxies within 0.5r„> was less than 0.5 (for 
the fossil groups this had to be larger than 2 magnitudes). The 
control sample comprises 3786 FoF galaxy halos. 

Figure [1] shows a plot of stellar masses versus FoF masses 
for both fossil (squares) and control (crosses) samples. To avoid 
infe rring results depende nt on the mass of the brightest galax- 
ies (iDe Lucia et al.ll2()06h from the control and the fossil sam- 
ples, we selected subsamples with approximately equal numbers 
of fossil and non-fossil groups, where the distribution of stellar 
masses of the central galaxies had been matched. This exercise 
provided two samples of fossil and non-fossil groups each with 
about 680 objects. We note that by insisting that the distributions 
of stellar masses of the central galaxies were identical for the 
two samples, the virial masses of the groups themselves were, 
on average, smaller for fossil than for non-fossil groups. This ef- 
fect can be seen in Fig. [T]if we match the y-axis distributions. 



However, this procedure is justified since we were interested in 
the properties of the first-ranked galaxies only. 

3. The brightest group galaxies 

3.1. The formation and evolution 

Given the availability of merger trees of galaxies in the 
Millennium Simulation, it is possible to s tudy the format i on and 
evo lution of central galaxies i n groups. iDe Lucia et al.l ( |2006|) 
and lDe Lucia & Blaizod (l2007h studied the evolution of galaxies 
of different stellar mass. They defined a set of particular times 
related to formation and evolution of galaxies. Following their 
work, we analyse the properties of the central galaxies of fossil 
and non-fossil groups, both samples having the same stellar mass 
distributions. Briefly, the different times are defined as follows: 

- Assembly time {za) is the time when 50% of the final stellar 
mass was already contained in a single galaxy. 

- Formation time (z/) is the time when half of the mass of 
the stars contained in the final galaxy at redshift zero have 
already formed. 

- Identity time (z,) is the time when the latest major merger 
occurred ( major merger is adopted as MilMi < 3 ( AI1 > 
M2)). 

- Extended identity time {ze) is the latest time when the sum of 
the masses of all progenitors merging at that time was greater 
than a third of the mass of the main progenitor (multiple si- 
multaneous minor mergers). 

We compute all of these times for the brightest galaxies in 
fossil and non-fossil groups, and also the assembly time of the 
FoF halos. In Table|2] we show the median of the different char- 
acteristic times and also include the probabilities that both dis- 
tributions, for fossil and non-fossil groups, are drawn from the 
same distribution (Kolmogorov-Smirnov test (K-S test)). K-S 
tests indicate significant differences in all cases. However, given 
the shape of the distributions, in some cases the median value is 
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Fig. 2. Cumulative distribution of assembly, formation, identity, and extended identity times (from top to bottom). Thick lines 
correspond to fossil groups, while thin lines are for non-fossil groups. 



Table 2. Median redshifts and probability values from K-S test: 
fossil versus non-fossils 





Fossil 


Non-Fossil 


K-S 


Za 


0.5642 


0.6871 


6 X IQ-" 


z/ 


3.5759 


3.5759 


6 X IQ-* 


Ze 


0.3197 


0.4566 


5 X IQ-* 


Zi 


0.2798 


0.4566 


1 X IQ-' 


ZaiFoF) 


0.4566 


0.1749 






insufficient to show these differences (for instance, when the me- 
dian of Zf in both samples is the same), therefore, in Fig.|2] we 
show the cumulative distributions of times in Fig.|2] Thick lines 
correspond to fossil groups, while thin lines are for non-fossil 
groups. Brightest galaxies of fossil groups have (a) assembled, 
(b) formed their stars, (c) experienced their last major merger 
and (d) experienced multiple simultaneous mergers, all at lower 
redshifts than bright central galaxies in non-fossil groups. In ad- 
dition, we confirm that the FoF halos of fossil groups assembled 
earlier than non-fossil groups. 

To check if the fact that we are analysing fossil and non- 
fossil groups with different halo mass ranges could introduce 
any bias, we repeated the process described above by matching 
the distributions of group virial masses for fossil and non-fossil 
groups. W e then succeeded in r eproducing all of the results de- 
scribed bv lDariush et alj ( |2007|) . in particular, that fossil groups 
in general, assembled earlier than non-fossil groups, for a given 



range of group virial masses. We also found that central galaxies 
in fossil groups have assembled, and experienced both their last 
major merger and multiple mergers all at lower redshifts than the 
central galaxies in non-fossil groups, but they have formed their 
stars at higher redshifts. This result resembles the analysis of 
iDe Lucia et alj ( |2006|) since by matching the virial masses of the 
groups, the resulting sample of fossil groups have central galax- 
ies with stellar masses typically larger than those of non-fossil 
groups. 

Regardless of whether either the stellar masses of the central 
galaxies or the virial masses of the halos are matched, although 
the FoF galaxy halos of fossil systeni s have assembled earlier 
dPariush et alj ( l2007[) : lD^nghia et all ( l2005l) . confirmed in this 
work using an independent semi-analytic model), their central 
galaxies have, on average, assembled later and, even more im- 
portantly, first-ranked galaxies of fossil groups also continue to 
experience major mergers for a longer period of time than in 
non-fossil systems. 

3.2. The morphologies and morphological mix during the last 
major merger 

The criteria used to select fossil groups did not include any con- 
straints on the morphological type of the first-ranked galaxy. 
However, it is interesting to study whether the first-ranked 
galaxy has any tendency to be an early-type galaxy. 

We must consider the shortcomings of the semi-analytic 
models when we attempt this analysis given that, historically. 
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most semi-analytic models have failed to reproduce some trends, 
although the new versions are fairly good at reproducin g many 
obser vational results related to morphology (see Bert one et al.l 
l2007h . However, it is still interesting to observe gene r al tren ds, 
defining morphology as completed bv lDe Lucia et alj (|2006|) . A 
galaxy is classified as elUptical if AM < 0.4 (AM = Mbuige - 
Miotai in the B-band), spiral if AM > 1.56 and SO in between. 
These authors state that this morphological type determination 
is robust for galaxies with stellar masses larger than a few times 
10^ h^^Aio- As can be observed from Table [1] we are dealing 
with galaxies beyond this limit. We find that 12% of central 
galaxies of fossil groups in the MSGC are non-elliptical (5% 
S and 7% SO's). Similar percentages are found in the control 
sample (3% are S and 9% are SO). 

Regarding the luminosities of the first-ranked galaxies in fos- 
sil groups, we found that even when no specific selection crite- 
rion was used to select luminous, central, group galaxies, the 
final sample of fossil groups had quite a bright median magni- 
tude of ~ -23.5. W e note that a search for fossil groups in the 
photometric SDSS (ISantos et al.ll2007l) . using similar selection 
criteria (but with no lower limit to the mass of the groups) pro- 
duced a sample of groups with a mean absolute magnitude of 
Mr -5 log h = -23.74 

iKhosroshahi et al.l (l2006l) found that central galaxies of fos- 
sil groups are different in their isophotal shapes (which are often 
disky) compared with the central galaxies of non-fossil systems 
(which often have boxy shapes). They then suggested that the 
central galaxies of fossil groups could be the result of wet merg- 
ers (gas-rich mergers) unlike the galaxies in non-fossil groups. It 
is, therefore, interesting to study the morphologies of objects that 
represent the pre-merging subclumps that produces the central 
galaxy, after a major merger. We call the process a wet merger 
if at least one of the two subclumps is a spiral galaxy, a dry 
merger if both subclumps are elliptical galaxies and mixed if 
one of them is elliptical and the other is an SO or both are SOs. 
Our results are quoted in Table[3] In contrast to expectations, we 
find that the vast majority of central galaxies in fossil groups in 
the MSGC are not produced by wet major mergers. In fact, cen- 
tral galaxies in both fossil and non-fossil groups appear to have 
undergone the same type of mergers, i.e. from the same morpho- 
logical mix, and are mostly the result of gas-poor mergers. This 
analysis was repeated for the subsample of systems that have el- 
lipticals at their centres (88% of the total sample) and similar 
results were found (val ues are also quoted in Table [3]l. 

On the other hand, iKhochfar & Burker3 (l2005h investigated 
whether the observational isophotal shape distribution of ellipti- 
cal galaxies could be reproduced in semi-analytic models by us- 
ing both the morphology of galaxies that merge in the last major 
merger and their mass ratio. Following their work, we classified 
the central elliptical galaxies that have had a major merger into 
"boxy" or "disky" according to the following criteria: 

- Ellipticals that experience last major mergers between two 
bulge-dominated galaxies (Albuige ^ 0.6 Altot) produce 
"boxy" remnants independently of the mass ratio 

- Last major mergers with mass ratio 1 < Mi /M2 < 2 produce 
"boxy" ellipticals. 

- Last major mergers with mass ratio 2 < Mi /M2 < 3 produce 
"disky" ellipticals. 

Our results for the sample with elliptical galaxies at their centres 
that have had a major merger are quoted in Table [3] This result 
disagrees with the predictions of Khosroshahi et al. (2006) based 
on 7 elliptical galaxies in fossil groups. We found no differences 
between the isophotal shapes of elliptical galaxies in fossil or 



Table 3. Morphological mix of subclumps that have merged to 
form the brightest central galaxy 





Full Sample 




Fossils 


non-Fossils 


DRY (E+E) 


48% 


39% 


Mixed (E+SO or SO+SO) 


22% 


26% 


WET (S+E or S+SO or S+S) 


8% 


9% 


No major merger 


22% 


26% 




E galaxy at the centre 




Fossils 


non-Fossils 


DRY 


57% 


45% 


Mixed 


26% 


31% 


WET 


8% 


9% 


No major merger 


9% 


15% 


boxy 


76% 


75% 


disky 


24% 


25% 



non-fossil systems, and in both the likelihood of being boxy was 
higher than being disky. However, the particular se mi-analytic 
model used in this work jPe Lucia & Blaizo3l2007[) could pro- 
duce a slightly higher fraction of boxy ellipticals than observa- 
tions (for instance, for non-fossil groups we found that for a 
median B-absolute magnitude of -21 , the ratio Nb o xvlNdi skv is 
around 2.3 while the observations of Bender et al.l d 1 9921) and 
the predictio ns of Khochfar & Burkert (|2005|) are close to 1.5 
(see Fig. 3 in lKhochfar & BurkertI ( l2005l) ). 



4. Mock galaxy catalogue and observations 

Although the principal aim of this paper is to study the 
brightest galaxies of fossil groups identified in the Millennium 
Simulation, a direct comparison of the fossil groups identified in 
the both simulations and observational catalogues can be used as 
a powerful test. We have particular interest to compare the frac- 
tion of FoF systems that are fossil groups in both a mock cat- 
alogue constructed from the MSGC and the Sloan Digital Sky 
Survey. 

To derive results that can be compared directly with obser- 
vations, we construct a mock catalogue in redshift space us- 
ing a snapshot at z = of the MSGC. We compute observer 
frame galaxy apparent magnitudes from the rest-frame absolute 
magnitudes provided by the semi-analytical model and tabu- 
lated (k+e)-corrections from iPoggiantil (Il997h . where the cor- 
rections were calculated according to an evolutionary synthe- 
sis model that reproduces the integrated galaxy spectrum in the 
range 1000 - 25000A. Distorted redshifts are computed using 
the pecuUar velocities of each galaxy. The mock catalogue com- 
prises ~ 1.3 X 10^ galaxies with apparent magnitudes lower than 
R - 17. 77 within the volume of one simulation box (zmo r ~ 0.17, 
7r/2 sr). 

Groups of galaxies in the mock catalogues (FoF galaxy ha- 
los) were i dentified by using an al gorithm similar to that de- 
veloped by iHuchra & Gellerl (Il982h . which is an adaptation of 
the FoF algorithm that takes account of the distortion caused 
by peculiar motions (redshift space) and the apparent magnitude 
cut-offs. As performed in real space, groups were identified as 
clusters with an overdensity of 200 and above. Particular care 
was taken in estimating the group centres. We used a method 
that computed the projected centre positions by weighting ap- 
propriately thepositions by the local densities and luminosities 
(iDiaz et al.ll2005h . 
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The masses, virial radii, and velocity dispersions of the FoF 
galaxy group were computed using Equation[Tl although the ra- 
dial velocity dispersions (cr,,) where then calculated using the bi- 
weight estimator described bv lBeers et al.lll990l (cr^n = V3cr,,). 
Only FoF galaxy halos with more than ten members were con- 
sidered. 

4.1. Fossil Groups in a mock catalogue 

The criteria used to select fossil groups was described in 
Sect. 12.31 Besides considering only groups with masses larger 
than 5 x 10'^ /i 'Al©, we also restricted the depth of the FoF 
galaxy halos by selecting groups with redshifts lower than 0.1. 
This restriction was based on statistical analyses that revealed 
that the reliability of the identification algorithm in redshift 
space increases for groups at redshifts below z/,„, =0.1. 

We found a comparatively low number of fossil groups: 22 
FoF galaxy halos were classified as fossils, which represented 
~ 3% of the FoF galaxy groups in the mock catalogue. Median 
properties of these groups were also quoted in Table [1] For the 
mock sample, we estimated the mean density of fossil groups 
up to the median redshift of the sample, which is a fairly robust 
measure of density. We measur ed 1.4x IQ-^h^ Mpc-\ 

It can be observed that the fraction of groups that satisfy the 
fossil criteria in the mock catalogue is lower than the ~ 5.5% 
found in MSGC. This discrepancy occurs for the following rea- 
son. In real space (MSGC) the search for galaxies whose mag- 
nitudes respect the AM12 > 2 selection criterion is completed 
within a sphere, which cannot be replicated for the mock cat- 
alogue given that we operate in redshift space, and by then 
measuring projected distances. When considering projected dis- 
tances in the selection of fossil groups, bonafide members in the 
outskirts of the groups may fail to satisfy the AM12 magnitude 
criterion, which then produces a smaller fraction of fossil groups 
(a similar effect is seen when fossil groups are defined to be 
within lr,,„. instead of 0.5rv,>). This hypothesis was tested in the 
mock catalogue by reidentifying fossil groups in the mock cata- 
logue, but, instead of using redshift space information, we used 
the available real space information to search for galaxies within 
a sphere of radius 0.5r,,„ . The percentages that we found in this 
simple test were far more similar to those drawn from MSGC. 

We note that in observational catalogues the entire identi- 
fication procedure will be affected in a similar way as that of 
our mock catalogue. We must therefore recall that percentages 
or number densities drawn from observations will be underes- 
timated. In particular, in our study case, we derived 3% from 
the mock catalogue (which should be equivalent to the obser- 
vations) but the real number should be closer to 5.5% for the 
fraction of groups that are fossils, from the entire population of 
groups/clusters with masses larger than 5 x lO'^ /z ' Al©- 

4.2. Comparison with SDSS 

We used the main galaxy sample with spectroscopic redshifts of 
the SDSS DR6 to identify observational fossil groups by using 
the same methodology applied in our mock catalogue. This sam- 
ple consisted of 574701 galaxies with apparent magnitudes in 
the r-band lower than 17.77 and a median redshift of 0. 1 . Briefly, 
FoF galax y halos were identified b y usi ng an algorithm similar 
to tha t of lHuchra& Gelled (Il982h (see iMerchan & Zandivare j 
I2OO5I for further details of the procedure) with an overdensity 
contrast of 200. Then, halos with more than 10 members, masses 
larger than 5 x 1 ' ^ /i" ' Me and redshifts lower than 0. 1 were se- 



Table 4. Percentage of Groups that are fossils in four different 
mock catalogues 





with close pairs 


without close pairs 


R < 17.77 


3% 


4% 


14.5 <R< 17.77 


0.2% 


0.26% 



Notes: Mock with close pairs: described in sect|4| Mock with- 
out close pairs: reproducing the fibre collision effect. Both were 
performed with different magnitude cut-offs. 



lected. Within this sample, we searched for fossil groups, i.e. 
groups with a magnitude gap of larger than 2 magnitudes in the 
r-band when considering galaxies within 0.5r,,„.. We found 6 fos- 
sil groups in this sample, which represented a fraction of 0.55% 
of all groups with more than 10 members and masses larger than 
5 X 10'^ /i^'AIq. Median properties of these fossil groups are 
quoted in Table [T] 

At first sight, this result is surprising, given that we predicted 
in the previous section that this fraction should be close to 3%. 
To understand this low percentage of fossil groups in the SDSS, 
we carefully considered the causes of incompleteness in the 
SDSS and we then reproduce these incompleteness in our mock 
catalogue. Two problems were quickly identified: 1) the incom- 
pleteness caused by fibre colUsions and 2) the incompleteness 
caused by the fibre magnitude limits and the image deblending 
software. The former meant that a fraction of ~ 70% of galaxies 
that have a neighbour closer than 55" were missing. The latter 
caused the spectroscopic sample to become noticeably incom- 
plete for galaxies brighter than r - 14.5. We then performed four 
different mock catalogues considering different combinations of 
these incompleteness, and we report in Table |4] the percentages 
of groups considered to be fossils. It can be seen that the fibre 
collision effect introduces a small bias in the percentages of fos- 
sil groups identified (which can be noted by comparing columns 
"with close pairs" and "without close pairs", in Table |4| but, in 
contrast, the lack of galaxies brighter than 14.5 strongly bias the 
results, causing the fractions to be far lower than the expected 
3-4% predicted from the full mock catalogue. Taking the incom- 
pleteness into account, the mean density of mock fossil groups 
up to the median redshift of the sample is 1 .2 x lO^^/i^ Mpc^^ 

To check this result, we repeated our search for groups in the 
SDSS, completed above, but now considering only galaxies in 
the magnitude range 14.5 < r < 17.77. We found exactly what 
we expected from the simulations: 0.27% of the FoF galaxy ha- 
los were fossil systems (this number can be directly compared 
with the final line in the final column of Table|4|. For the remain- 
ing fossil groups, the mean density up to the median redshift was 
1.5 X 10^^ Mpc^^, which was quite similar to that obtained in 
the mock catalogue after all biases were considered. 

The lack of bright galaxies also caused the few fossil groups 
that we identified in the SDSS DR6 to have first-ranked galax- 
ies that were not as bright as those found in the semi-analytical 
model. 

Despite the small number of fossil groups identified in the 
main spectroscopic galaxy sample of SDSS, we attempted to 
analyse their main properties. In Table |5] the properties of 
each of the six groups are quoted. Figure [3] shows the angu- 
lar distribution of the equatorial coordinates a versus 6 for all 
known fossil groups. Grey points are the galaxies in the main 
sample of the SDSS. Crosses are the six fossil groups pre- 
sented in this work. Filled circles are the fossil groups listed in 



Table 5. Fossil groups in the main galaxy sample of the SDSS DR6 



name 


[h:m:s] 


Sen 

[d:m:s] 






O-iD 

[km .r'] 


t'vir 


[h:m:s] 


<5i 
[d:m:s] 




Ml - 5log{h) 
[r-band] 


eClassi 


I 


07:35:29.04 


39:45:21.6 


0.0866 


11.3 


797 


0.86 


07:35:36.46 


39:45:53.57 


0.0874 


-22.43 


-0.190478 


II 


07:58:34.32 


37:46:12.0 


0.0406 


7.5 


563 


1.01 


07:58:28.11 


37:47:11.87 


0.0408 


-22.29 


-0.151681 


III 


15:19:02.64 


04:17:49.2 


0.0466 


8.1 


579 


1.04 


15:19:03.52 


04:20:01.14 


0.0468 


-22.09 


-0.177914 


IV 


07:34:33.84 


26:53:56.4 


0.0796 


7.9 


551 


1.12 


07:34:22.22 


26:51:44.93 


0.0797 


-22.35 


-0.197101 


V 


14:38:55.68 


10:08:16.8 


0.0552 


1.1 


523 


1.21 


14:38:47.60 


10:07:17.26 


0.0553 


-21.70 


-0.158786 


VI** 


08:56:37.44 


21:52:51.6 


0.0823 


9.9 


544 


1.43 


08:56:51.80 


21:49:49.44 


0.0827 


-22.17 


-0.175508 



Notes: q'„„,<5„„: equatorial coordinates of the centre of mass of the group; Zc,,,'- spectroscopic redshift of the centre of mass; Mvi,-'- virial mass of the group; 
ctsd: 3-D velocity dispersion of the group; r,.,,-: virial radius of the group; a\. Si, zi equa torial coordinates an d spectroscopic redshift of the brightest galaxy; 
Ml — 5log{h): rest frame absolute magnitude of the brightest galaxy, k-l-e corrections from lBlanton et al.l2003l ; eClassi : spectral type parameter of the brightest 
galaxy, provided by SDSS. 
** It is not a fossil group (see text) 
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Fig. 3. Grey points are the galaxies of the m ain sample of the SPSS. Cross es are the six fossil groups presented in this work. 
Filled circle s are the fossil grou ps found in the lMendes de Oliveira et ani2006l with z < 0.1, and empty circles are the fossil groups 
identified bv lSantos et"ani2007l with z < 0.1 



iMendes de Oliveira et al ] (l2006h with z < 0.1 and em pty circles 
are the fossil groups identified bv lSantos et al.l(l2007l) in the pho- 
tometric sample of SDSS with z < 0.1. There are three filled cir- 
cles in the SDSS area that we failed to identify. Using the SDSS 
DR6 Finding Chart Tool0, we analysed a field about the centre 
of those fossil groups and realised that they could not be identi- 
fied because of different reasons. RX J 133 1.5+1 108 could not be 
identified because its main galaxy could not be observed due to 
possible fibre collisions. We failed to identify RX Jl 159.8+5531 
as a fossil group because the two-magnitude gap was not present 
in this group. We identified this group instead to be massive 
and containing numerous galaxies, but the two-magnitude cri- 
terion caused us to reject this halo (the absolute magnitude dif- 
ference in R-band between the first and second ranked galaxy is 
AMi2 = 1.6). Finally, NGC 6034 was not detected since we did 
not identify a massive and numerous halo in its position. Instead, 
we identified the central galaxy to be part of the outskirts of a 
larger group. A closer inspection also indicated that this galaxy 
also had a close bright neighbour in the same large group, caus- 
ing the magnitude gap to be AM12 - 1 .07 and not 2, as the strict 
definition of fossil groups requires. 

It is also int e restin g to compare the fossil groups identified 
by ISantos et al.l (|2007|) and ourselves. Only one of six groups 
(number II in our Table |5]l was in their catalogue. One reason 
for these differences is that lSantos et al.l (|2007|) insisted that all 
fossil groups had extended X-ray emission in the ROS AT all-sky 



catalogue. For our six groups, one has extended X-ray emission 
(group II) and one has point source X-ray emission (as detected 
by XMM, group number IV) and four do not have any detectable 
ROSAT emission (and have not been observed by other satel- 
lites). Given that the sensitivity of ROSAT is quite low, these 
groups should be studied again with deeper X-ray observations 
to confirm their nature as fossil groups. Another important dis- 
crepancy between the two works is the way in which a fossil- 
system is defined. Differences in the radii used to test the mag- 
nitude gap criterion (virial radius or fixed radius in kpc) and the 
way of defining the group membership may cause each study to 
select different systems. 

By using the NASA/IPAC Extragalactic Database 0, we also 
analysed whether the SDSS fossil groups presented in this work 
had been previously identified as groups or clusters of galaxies. 
Fossils I, IV, V, and VI did not have previous identifications as 
systems of galaxies. Fossil II correspond ed to the X-ray cluster 
NGC 2484 GROUP ( iPopesso et alj|200l and fossil III was pre - 
viously identified as the cluster CAN 245 dWegner et al.lll999l) . 

We also tested the possibility that our fossil groups were spu- 
rious identifications due to the incompleteness caused by the fi- 
bre collisions or the fibre magnitude limit that affect the spec- 
troscopic sample of the SDSS. Using the SDSS DR6 Finding 
Chart Tool, we analysed a field of 0.5 r,.,> around the centre of 
each fossil group and found that in fossil VI the brightest galaxy 
of this group has a close neighbour that is a spectroscopic tar- 



^ http://cas.sdss. org/dr6/en/tool s/chart/chart . asp 



^ http://nedwww.ipac.caltech.edu 
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get not observed due to the fibre collision problem. The mag- 
nitude gap between these two galaxies is less than two magni- 
tudes (AM12 ~ 1 .33). Although the redshift of the second-ranked 
galaxy is unknown, we note that it is likely that this is not a fossil 
group. This is one of the groups that are not detected in X-rays. 
A similar situation applies to fossil group III, since in the area 
corresponding to half the virial radius there is a galaxy whose 
(apparent) magnitude differs by less than 2 from the (apparent) 
magnitude of the brightest galaxy, but this galaxy was missing 
from the SDSS survey because it has a very close neighbour. 
However, in this case, the visual inspection is insufficient to con- 
firm whether this galaxy belongs to the same group. 

In summary, the spectroscopic sample of SDSS DR6 is 
,probably, not the most appropriate for the of study observational 
fossil groups due to its inherent incompleteness, which means 
that we are often unable to or falsely identify groups. The enor- 
mity of this galaxy data set is still however a useful guide for 
deeper optical observational studies or even X-ray observations, 
which are beyond the scope of this paper 

5. Conclusions 

To help interpret observations, we have analysed the properties 
of the first-ranked galaxies of simulated fossil groups and their 
main differences with respect to the brightest galaxies of non- 
fossil systems. To perform this analysis, we have studied a sam- 
ple of fossil groups obtained from the largest simulated galaxy 
catalogue at present, the Millennium simul ation combined with 
a sem i-analytic model of galaxy formation (iDe Lucia & Blaizoj 
I2OO7I) . The sample of fossil groups identified was sufficiently 
large to enable robust statistical analyses to be completed. A 
control sample of non-fossil groups was also selected. The goal 
of this research was twofold: to predict results based on semi- 
analytical models about fossil groups and test the semi-analytical 
model of formation and evolution by comparing with observa- 
tional results. 

Our analysis was performed in both the Millennium Galaxy 
Catalogue (MSGC) and a mock catalogue. The fraction of fossils 
predicted by the mock catalogue was compared with the fraction 
of fossils found in the SDSS data set by applying the same al- 
gorithm of identification. If the two main sources of incomplete- 
ness affecting the spectroscopic sample of SDSS were taken into 
account, the fraction of fossil groups found in the mock cata- 
logue was consistent with that observed in the SDSS. We have 
presented six candidate fossil systems in SDSS DR6. One has 
been rejected after visually analysing the surroundings of the 
brightest galaxy. The remaining groups await X-rays observa- 
tions or deeper photometric studies to confirm their nature as 
fossil groups. 

Our main results can be summarised as follows: we confirm 
the old age feature of fossil systems and that fossil groups with 
masses larger than 5 x lO'-' /j^'AIq represent ~ 5.5% of mas- 
sive systems in the same mass range, in the Millennium Galaxy 
Catalogue and a lower percentage (of 3%) of similar systems in 
a mock galaxy catalogue. Fossils in the Mock Catalogue were 
identified in redshift space (just as achieved in observations), 
where galaxies from the outskirts of the groups affect the AM12 
selection criterion and reduce the resulting number of fossils. We 
found that 88% of fossil systems have central galaxies that are 
ellipticals. In addition, we found that the first-ranked galaxies of 
fossil groups and non-fossil groups have the same morpholog- 
ical mixtures, which can be considered to be the results of gas 
poor mergers. If the central galaxy of a fossil group is always el- 
liptical, the likelihood that its progenitor galaxies merge in a dry 



merger is higher Although this result disag rees with the obser- 
vational studies of iKhosroshahi et all (l2006h for seven elliptical 
galaxies in fossil groups, further observational data is required to 
fully resolve this issue. Finally, we investigated the nature of the 
central galaxies by keeping track of their merging histories. On 
the one hand, the fossil groups in general assembled earlier than 
non-fossil groups, whilst their central galaxies assembled later 
We also found that first-ranked galaxies in fossil groups have 
undergone a major merger later than their counterparts in non- 
fossil systems. We expect this result to be confirmed by future 
observational catalogues. 
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